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ABSTRACT 


This  paper  presents  a  "one  pass"  algorithm  that  determines  an  "advanced" 
dual  basic  feasible  solution  for  a  class  of  capacitated  generalized  network 
problems.  Special  cases  in  this  class  of  problems  include  transportation  and 
transshipment  problems.  Computational  results  are  included  which  show  that 
this  new  start  substantially  improves  the  solution  performance  of  the  dual 
method  for  transportation  and  transshipment  problems.  In  fact,  a  dual  code 
employing  this  advanced  start  is  found  to  be  faster  (in  terms  of  total  solution 
time)  than  the  fastest  out -of- kilter  code  SUPERK  on  highly  rectangular  trans¬ 
portation  problems. 


II / 


Introduction 
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This  paper  presents  a  "one  pass"  algorithm  that  determines  an  "ad¬ 
vanced"  dual  basic  feasible  solution  for  a  class  of  capacitated  generalized 
network  problems.  Special  cases  in  this  class  of  problems  include  the 
standard  transportation  problem  and  the  pure  network  or  transshipment  prob¬ 
lem. 

The  algorithm  presented  is  "advanced"  in  that  it  determines  a  solution 
which  is  not  onlj  feasible  but  relatively  close  to  the  optimum  of  the  dual 
of  a  generalized  network  problem.  The  advanced  algorithm  is  a  direct  ex¬ 
tension  of  the  dual  start  algorithm  of  Glover,  Klingman  and  Napier  [ 13  3  - 
Their  algorithm  was  the  first  to  obtain  dual  feasible  solutions  for  a  class 
of  multi-source,  multi-sink,  capacitated  generalized  network  problems,  and 
extended  the  one  pass  algorithm  of  Charnes  and  Raike  [4j  for  determining 
an  optimal  solution  for  another  class  of  generalized  network  problems  The 
advanced  algorithm  of  this  paper  differs  from  the  earlier  dual  start  [13] 
in  that  objective  function  information  is  used  in  constructing  the  dual 
feasible  solution. 

This  research  was  motivated  by  recent  computational  investigations 
[1,  9,  11,  12,  17,  19]  which  have  found  that  special  dual  simplex  approaches 
to  transshipment  and  transportation  problems  are  generally  inferior  to 
special  purpose  primal  simplex  and  primal  dual  (or  out-of-kilter)  approaches 
on  large  problems.  In  particular  these  studies  have  found  that  the  special 
purpose  primal  simplex  approaches  [9>H>17]  are  substantially  (^-10  times) 
faster  than  out-of-kilter  approaches  and  that  the  out-of-kilter  codes  are 
faster  than  dual  codes  which  employ  the  dual  start  in  [133  *  Since  the 
study  [1]  showed  that  it  is  possible  to  increase  the  efficiency  of  the  out- 
of-kilter  method  by  at  least  a  factor  of  four,  we  felt  that  further  com¬ 
putational  investigation  on  dual  approaches  was  warranted.  Furthermore, 
since  the  reported  computational  testing  on  the  dual  approach  [9>  12] 
indicated  that  its  major  computational  weakness  lies  in  the  large  time  re¬ 
quired  to  execute  a  pivot,  it  seemed  reasonable  that  starting  procedures 
which  would  reduce  the  number  of  pivots  required  to  find  an  optimal  solution 
should  be  designed  and  computationally  tested. 

Computational  investigation  on  the  advanced  start  developed  herein 
shows  that  it  acheives  this  reduction.  In  particular,  the  results  indicate 
that  the  advanced  start  reduces  both  the  solution  time  and  number  of  pivots 
over  other  auaj.  upproeches  using  the  start  in  [133  hy  roughly 
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While  this  result  still  is  not  able  to  make  the  dual  approach  competitive 
with  the  primal  approach  on  most  transportation  and  transshipment  prob¬ 
lems,  our  computational  testing  does  indicate  that  for  highly  rectangular 
transportation  problems  (i.e.,  transportation  problems  with  only  a  few 
origins  and  a  large  number  of  destinations)  the  dual  approach  using  the 
advanced  start  procedure  is  superior  to  the  fastest  known  out-of-kilter 
code  [1]  and  competitive  with  the  fastest  primal  codes  [9>  11,  17]. 

A  principal  motivation  cited  by  other  papers  [3>  20]  lor  applying  the 
dual  method  to  network  problems  is  the  possibility  that  the  problem's 
'supplies'  and  demands'  may  not  be  permanently  fixed,  but  rather  are  sub¬ 
ject  to  change.  In  such  a  situation,  the  ability  to  begin  from  an  optimal 
bas  s  to  a  given  problem  and  proceed  via  the  dual  method  to  an  optimal  sol¬ 
ution  for  the  altered  problem  is  extremely  useful.  Our  computational  results 
indicate  that  the  advanced  start  procedure,  however,  offers  a  new  motivation 
f  _>r  the  use  of  the  dual  method — namely,  by  use  of  our  algorithm,  it  ic  pos¬ 
sible  to  employ  the  dual  method  not  only  for  'postoptimization, '  but  also 
for  solving  the  original  problem  directly. 


2.  Problem  Statement 

A  transshipment  problem  will  be  defined  as  consisting  of  m  nodes  that 
are  connected  pairwise  by  a  collection  of  n  directed  arcs.  Let  N  denote 
the  set  of  all  admissible  arcs  (i,j)  in  a  transshipment  problem,  then  a 
capacitated  generalized  transshipment  problem  can  be  stated  mathematically 
as: 

minimize  2  c. .  x  . . 

(i,j)€N  1J 

subject  to  -2  x  +  2  k..x.  *  b.,  (i  =  1,2, ...,m)  (l) 

1}  3  1 

Wtti4  ’  [(i’j)cNl 

where 

1.  x. .  is  the  flow  from  node  i  to  node  j. 

ij 

2.  c„  is  the  cost  of  sending  a  single  unit  of  flow  from  node  i  to  node  j. 

3.  b^  is  the  amount  of  supply  (demand)  at  node  i,  where  the  supply  (demand) 
at  node  i  is  denoted  by  a  negative  (positive)  b^ 

h.  k  is  the  attenuation  factor  affecting  the  amount  of  flow  leaving  the 

ij 

initial  node  i  and  flowing  into  terminal  node  j. 
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The  dueil  to  Problem  (l)  may  be  stated  as: 

m 


maximize 


£  b  w.  +  L  s.  L 

i*l  1  1  (i*j)eN  J  1J 


L 


subject  to 


k  w .  -  w.  +  s .  . 
ij  i  1 


t.  .  =  c,  . 
ij  io 


unrestricted 

s.  .  t  .  >  0 
ij  iJ  ~ 


Ui,i)  €  N]  (2) 

(i  =  1,2, ... ,m) 
e  N] 


The  advanced  algorithm  presented  in  the  next  section  finds  a  basic  feasible 
solution  (if  it  exists)  for  Problem  (2)  provided  (a)  c  >0,  (i,j)  c  N, 


(b)  0  <  k.  .  <  1, 
ij  - 


ij 

(i,j)  c  N,  and  (c)  the  rank  of  the  incidence  matrix 


is  in  -  1. 


Conditions  (a)  -  (c)  are  not  as  restrictive  as  they  may  appear. 
Requirement  (a)  is  generally  innocuous  in  operations  research  applications 
[22,  p.  178].  Restriction  (b)  prohibits  the  inclusion  of  a  loop  in 
the  dual  feasible  basis  determined  by  the  algorithm.  Condition  (c) 
is  always  satisfied  by  transportation  ar.d  pure  network  problems;  it 
is  only  restrictive  in  the  generalized  network  case  when  the  rank  of 
the  incidence  mt  trix  ism.  Even  in  this  case,  the  algorithm  will  generate 
a  dual  feasible  spanning  tree  solution,  although  not  a  basic  solution. 
Requirement  (c)  can  be  circumvented  if  the  generalized  network  contains 
a  single  master  source  With  an  unlimited  supply. 

In  order  to  genert.te  a  dual  basic  feasible  solution  to  Problem  (2), 
the  algorithm  first  derives  a,  basic  feasible  solution  to  the  problem: 


maximize 


'm 

£ 

i=l 


w. 

1 


subject  to  k_  w^  -  w±  <  c_  [(i,j)  e  N]  (3) 

w^  unrestricted 

Problem  (3)  is  the  dual  to  Problem  (l)  in  uncapacitated  form. 

The  algorithm  presented  in  the  next  section  will  generate  a  basic 
feasible  solution  to  Proolem  (3).  Clearly  such  a  solution  will  yield  a 
feasible  solution  to  Problem  (2).  The  algorithm  proceeds  by  first 
selecting  an  appropriate  starting  node  and  determining  the  generalized 


longest  path  from  all  other  nodes  that  can  reach  the  starting  node  by 
a  directed  arc.  One  of  the  above  nodes  is  then  selected  (and  its  asso¬ 
ciated  due.l  variable  value  fixed)  in  such  a  manner  that  dual  feasibility 
is  preserved  and  the  objective  function  (at  this  iteration)  is  mcocimized. 
At  each  subsequent  stage  generalized  snortest  and  longest  paths  are 
determined  from  fixed  to  non-fixed  nodes;  one  node  is  again  selected 
whose  dual  variable  value  preserves  feasibility  and  yields  the  best 
return  for  the  objective  function.  The  process  is  reprated  until  all 
m  nodes  of  the  network  have  been  selected.  The  algorithm  is  a  "one-pass" 
method  in  the  sense  of  Charnes  and  Raike  [i+]  in  that  at  each  interation 
one  new  node  is  selected;  the  process  then  yields  a  succession  of  m  single 
steps  which  determine  the  solution. 

Algorithm 

Let 

T(p)  =  the  set  of  "fixed"  nodes  at  iteration  p. 

S  =  the  set  of  arcs  from  fixed  nodes  to  nonfixed  nodes. 

R  =  the  set  of  arcs  from  nonfixed  nodes  to  fixed  nodes 

B  =  the  set  of  arcs  currently  selected  to  be  in  the  basis 

1.  Select  the  starting  node  t  to  be  a  source  node.  Set  the  dual  vari¬ 
able,  w^.,  of  this  ne_e  equal  M  =  ^  c 


(i>j)  eN 

(For  the  pure  network  problem,  any  node  can  be  selected  arbitrarily 
as  the  starting  node  t  and  M  can  be  set  equal  o  zero. )  Initialize 
by  setting  =  ^  &ni  R  m  ^  Set  q  =  t;  in  general  q  denotes 

the  node  whose  dual  variable  was  last  set.  Set  p  =  1;  in  general  p 
denotes  the  current  iteration  of  the  algorithm.  Set  w^  (0)  =  <*>  for 
i  {  T(l)  and  w^  (0)  =  0  for  i  /  T(l). 


2.  If  S  =  { ( i , j ) c  N:  ie  T(p),  j/  T(p))  is  empty,  proceed  to  step  3; 
otherwise  for  each  node  j/(  T(p),  let 
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w1  (p)  = 


r 

min  ( w"T  (p-l). 


w  +  c  . 


)  if  (q>j)e  N 
if  (q,j)/  N 


3.  If  R=  ((i,j)r  N:  i/  T(p),  jr  T(p) }  is  empty,  proceed  to  step 
otherwise  for  each  node  jf  T(p),  let 


0  (max  (w.(p-l),  k.  w  -c.  )  if  (j,q)c  N 
w2(p)  J  ™ 

J  /  o . 


w  (p-l) 

.  J 


if  (i,q)/  N 


4.  If  SWR  is  empty  and  there  exists  a  node  ^t  T(p),  then  stop  for  the 
network  is  disconnected;  otherwise  if  S0  0  set  w^,=  min  (w^(p) } 
(where  f  equals  the  T(p)  yielding  the  minimum)  afad^ir  0  set 


0  set 


2 

w^= .^ax^ (p) )  (where  r  equals  the  j/  T(p)  yielding  the  maximum). 

Next  compute  the  "effective  net  stock"for  node  f  (if  S t  0)  as 

bi  =  b_  +  2  b  ./k_ .  +  2  k  b . 

f  f  (f, j)e  N  J  fj  (j,f)c  N  Jf  J 
it  T(p)  it  T(p) 


and  if  R/  0  compute 

b1  =  b  +  2  b./k  .  +  2  k.  b. 

r  r  (r,j)e  N  0  rj  (j,r)e  N  jr  J 
it  T(p)  it  T(p) 


Then  if  w^,b^  >  w^b^  ,  set  q  =  f  and  augment  T(p)  by  q;  otherwise 


set  q  =  r  and  augment  T(p)  by  q.  If  q  =  f  augment  B  Dy  the  arc 

(p,f)  for  which  w^,  =  wp+  cpf  .  otherwise  augment  B  by  the  arc 
1  k  * 

V 

(r,p)  for  which  w  =  k  w  -c  .  In  either  case,  if  there  is  more 
r  rp  p  rp 

than  one  arc  then  choose  any  such  arc.  Finally,  if  T(p)  contains 
all  the  nodes  the  algoritlim  is  completed;  otherwise  set  T(p+l)=  T(p), 
p  =  p+1,  and  return  to  step  2. 


1.  Validity 


In  order  to  establish  the  validity  of  the  algorithm,  it  is  necessary 
to  first  prove  the  following  lemma: 
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Lemma:  At  each  iteration  p, 

rain{w^(p) }  >  max{wi(p)) 
i/T(p)  i/  T(p) 

(i.e.,  if  w  and  w  are  both  defined  at  iteration  p,  then  w^,  >  w  ). 

f  r  x  —  r 


proof : 

The  result  is  established  by  induction.  By  the  lemma  in  [13],  we 

1  2 

can  assume  w.l.o.g.  that  all  w^(p)  and  w^(p)  at  every  stage  p  are  non¬ 
negative  if  the  first  dual  variable  is  set  equal  to  M.  Additionally, 
we  assume  that  the  nodes  have  ueen  renumbered  in  theorder  in  which 

their  w.'s  are  fixed. 

1 

(i)  For  the  first  iteration  (p=l),  we  may  assume  w.l.o.g.  that 

M  +  c 

w  and  w  are  defined  and  Lave  the  values  w  =  If  and 

w  =  k  ,M  -  c  ,  Since  c. .  >  0  for  all  (i,j)e  N  and  0  <  k. .  <  1 
r  rl  rl  ij  -  v  ij  - 

for  all  (i,j)e  N,  w„  >  w  and  thus 

f  ~  r 

1  2 
min{wi(p) }  >  max{wi(p)}. 

i/  T(l)  i/  T(l) 

(ii)  Now  assume  the  hypothesis  is  true  for  the  p  stage. 

Thus 

1  2 
min{wi(p)}  >  wp  >  max{wi  (p)}  (4) 

i j  T(p)  U  ?(p) 

and  for  each  i^  T(p) 


(min{w^(p), 

(w,(p) 


w 

_E 


k  . 
Pi 


if  (p,i)e  N 

if  (p>i)/  N 


and 


2/  «nax{w.(p),  k.  w  -  c.  } 

w.(p+l)  =•<  iv^/7  ip  p  ip 


Wi  (P) 


We  may  assume  w.l.o.g.  that  for  each  i/  T(p) 
wj  (p+1)  =  minfwj(p),  Wp  +  Cpi) 


and 


Pi 


2  2 

w.  (p+l)  =  max(w.(p),  k.  w  -  c.  ) 

i  ip  p  ip 


if  (i,p)c  N 
if  (i,p)/  N 
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Therefore, 

W  *r  C 

min  {w*(p+l)}~  min  {  w^(p),  -2-£ — 21  } 

H  T(P+1)  i/  T(p+1)  pi 

and 

2  •*  g 

max(w. (p+l) ) =  raax{w.(p),  k.  w  -c.  ). 
i'  iv^/7  ip  p  ip 

i /  T(p+1)  i f  T(p+1) 

To  establish  that  min  (w^  (p  +  l))  >  max  {  w^  (p+l)), 
if  T(p+i)  ic  T(p+i) 

it  suffices  to  show  that  the  following  cases  are  true: 

(a)  min  {wi(p))>max  (w^p)} 

if  T(p+1)  if  T(p+1) 

W  “f*  0 

(b)  min  (  _ £i  )  >  max  {kipwp  -  cip) 

k  • 

if  T(p+l)pi  if  T(p+l) 

(c)  min  (w^(p))  >  max  ^kipwp  "  cip) 

i f  T(p+1)  if.  T(p+1) 

(d)  min  (Wp+  Cpi)  >  max.  (w?(p))  . 

if -T(p+ljpi  if  T(p+1) 

Case  (a)  follows  directly  from  the  induction  hypothesis  since 

2  2 

max  (wjL(p) }  <  max  {wi  (p)) 

if  T(p+1)  if  T(p) 

and 

2  2 
min  (wi(p))  >  min  (p)). 

if  T(p+1)  i  f  T(p) 

Since,  by  the  lemma  in  [13],  w  >0,  and  c.  .  >  0  and  0  <  k.  .  <  1  for  (i,j)e  N, 

p  —  ~  ij  — 

case  (b)  is  clearly  true.  Cases  (c)  and  (d)  follow  from  expression  (k), 
c. .  >  0  for  (i,j)e  N,  and  0  <  k.  .  <  1  for  (i,j)e  N.  Therefore,  the  lemma 

follows  by  induction. 


Theorem:  The  algorithm  obtains  a  feasible  basis  for  Problem  (2)  or  determines 


that  the  network  associated  with  Problem  (2)  is  disconnected. 
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Proof:  The  algorithm  terminates  alter  at  most  m  iterations,  because  at  each 
iteration  T(p)  is  augmented  by  one  node  or  else  the  algorithm  stops  with  the 
assertion  that  a  feasible  basis  has  been  obtained  or  that  the  network  is  dis¬ 
connected. 

The  assertion  that  the  network  is  disconnected  occurs  when  S  and 
R  are  empty  and  there  exists  some  node  j/  T(p).  Under  these  conditions  the 
network  is  split  into  at  least  two  disconnected  parts. 

Ruling  out  disconnectedness,  it  must  be  demonstrated  that  the 

values  assigned  to  the  w.  are  feasible,  i.e.,  k. .  w.-  w.  <  c. .  for  all 

i  ij  o  l  -  ij 

(i,j)e  N.  Again  the  method  of  proof  is  induction.  Assume  the  nodes  are 

numbered  in  the  order  in  which  they  are  fixed.  At  iteration  r,  the  number 

of  nodes  in  T(r)  is  r.  Suppose  that  the  w.  values  for  these  nodes  represent 

a  feasible  solution  to  the  dual  subproblem  obtained  by  considering  all  nodes 

in  T(r)  and  all  links  joining  the  nodes  in  T(r).  We  shall  show  that  the 

process  of  augmenting  the  (r+l)st  node  to  T(r)  assigns  to  its  dual  variable 

w^+1  a  value  that  is  feasible  for  the  dual  subproblem  obtained  by  considering 

the  (r+l)  nodes  in  T(r+l),  thereby  yeilding  the  desired  inductive  result. 


(i)  To  begin  ,  let  r  =  1.  The  assignment  w^=  M  clearly  provides  a  feasible 
solution  to  the  subproblem  consisting  of  the  first  node.  For  r  =  2,  the 
possibilities  are  that  w^=  11  +  c12  or  wp=  k^  M  -  cpi . 

k12 

If  w^=  M  +  c12  ,  then  w^  >  w,  since  0  <  k.  .  <  1  and  c.  .  >  0 

2  -r -  2-1  ij  -  ij  - 

*12 

for  (i,j)e  N.  Thus,  the  values  of  w^  and  w^  satisfy  the  feasibil¬ 
ity  c  editions  not  only  for  the  arc  (1,2),  but  also  for  the  arc  (2,1)  if  it 
exists.  On  the  other  hand,  if  wp=  k^  M  -  cp^  then  w^  <  w^.  Therefore, 
the  values  of  w^  end  w^  satisfy  the  feasibility  conditions  not  only  for  the 
arc  (2,1)  but  also  for  the  arc  (1,2)  if  it  exists. 

(ii)  Now  assume  that  the  r  nodes  in  T(r)  have  node  values  which  yield  a 
dual  feasible  solution.  It  will  be  shown  that  adding  the  (r+l)8^  node  then 
yields  a  dual  feasible  solution. 

Let  p(r+l)  =  {j  ;  j  e  T(r),  (j,  r+l)  e  N)  and 

let  q(r+l)  =  (j  ;  j  e  T(r),  (r+l,  j)  e  N). 
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Theri  automatically  w  ,  satisfies  the  condition  k.  ,  w  , -  w .  <  e. 

r+1  j>r+l  r+1  j  -  j,r+l 

for  jc  p(r+l).  By  the  previous  lemma,  w  ,  >  max  {k  ,  .w.-  c  ,  .•> 

0  *  r+1  -  r+1,  j  o  r+1, j} 

j<-  q(r+l) 

and  therefore  the  condition  k  .w.-  w  .  <  c  for  je  q(r+l)  is  satisfied. 

j  r+x 


Likewise,  suppose  node  r+1  is  being  added  in  the  "reverse  direction" 

(i.e.  w  .=  max  {k  ,,  .w.-  c  ,,  .)).  Then  the  condition  k  .w.-  w  , 
r  1  Jf  q(r+l)rl,j  0  r+1>J  r+1>0  J  r+1 

<  c  ,  .  for  je  q(r+l)  is  automatically  satisfied.  Again  by  the  lemma 
—  r+f>J 

w  <  min  {wj+  C j , r  +i  },  the  condition  k_.  _  .,w  -  w,  <  c .  for 

r+1  >  p(r+ir  kj,rtl_  0  “  J,r+1 

jc  p(r+l)  is  satisfied.  This  completes  the  feasibility  proof. 


The  final  step  of  the  proof  is  to  show  that  the  feasible  solution 
thus  determined  by  the  algorithm  is  in  fact  basic  feasible,  i.e.,  the  arcs 
in  B  form  a  basis  when  the  algorithm  terminates.  First,  upon  termination  b 
contains  rn-1  arcs  because  one  new  arc  is  added  to  B  each  time  a  node  is  added 
to  T(p),  except  at  initialization.  The  arcs  in  B  clearly  span  the  network. 
Further,  the  arcs  in  B  contain  no  cycles.  To  verify  this,  assume  the  con¬ 
trary  and  let  (r,s)  be  the  first  arc  added  to  B  that  creates  a  cycle  with 
the  previous  arcs.  Then  there  must  be  an  arc  (r,j)e  B  and  an  arc  (i,s)c  B. 
But  this  is  impossible  when  (r, s)  is  added  to  B,  since  all  arcs  (i,j)e  B 
satisfy  ie  T(p)  and  jc  T(p),  whereas  to  augment  B  with  the  arc  (r,s)  either 
r  /  T(p)  or  s /  T(p).  This  completes  the  proof. 
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Computational  Results 

To  test  the  advanced  start  procedure,  the  dual  codes  developed  in  [9, 
12]  were  modified  to  use  this  start.  These  codes  use  the  double  pricing  pro¬ 
cedure  [12,1A]  for  determining  the  unique  updated  linear  equation  which 
expresses  the  basic  variable,  leaving  the  basis  in  a  pivot,  as  a  linear 
combination  of  the  current  nonbasic  variables  (for  a  complete  description 
see  [12,1k]).  The  important  aspect  is  that  the  double  pricing  procedure  is 
quite  efficient  for  calculating  these  coefficients.  Thus  these  codes  should 
be  efficient  in  the  pivot  phase  of  the  simplex  dual  method  [20]. 

All  of  the  codes  used  in  this  study  are  in-core  codes;  i.e.,  the 
program  and  all  of  the  problem  data  simultaneously  reside  in  fast-access  mem¬ 
ory.  They  are  all  codec  in  FORTRAN  and  none  of  them  have  been  tuned  for  a 
par+icular  compiler.  All  of  the  problems  were  solved  on  the  CDC  6600  at  the 
Ur-.versity  of  Texas  Computation  Center  using  the  Run  compiler.  The  computer 
jobs  were  executed  during  periods  when  the  machine  load  was  approximately 
the  same,  and  all  solution  times  are  exclusive  of  input  and  output,  i.e., 
the  total  time  spent  solving  the  problem  was  recorded  by  calling  a  Real 
Time  Clock  upon  starting  to  solve  the  problem  and  again  when  the  solution 
was  obtained. 

Initially,  the  advanced  dual  start  code  was  tested  on  eighteen  of 
the  benchmarked  problems  in  [19]  since  these  problems  had  already  been 
solved  using  the  codes  by  Benninyton  [2],  Boeing,  General  Motors,  SHARE 
[5,21],  SUPERK  [l],  Texas  Water  Development  Board,  PNET-I  [9],  and  DNET 
[9],  The  fiist  five  of  these  codes  are  variants  of  the  out-of-kilter  method 
[7,8]  except  for  Bennington  [2].  PNET-I  [9]  is  the  fastest  known  special 
purpose  primal  simplex  code  for  solving  capacitated  and  uncapacitated  trans¬ 
shipment  problems.  DNET  is  a  special  purpose  dual  simplex  code  using  the 
start  procedure  in  [13]  for  solving  transshipment  problems.  DNET  was  modi¬ 
fied  to  use  the  advanced  dual  start  and  this  version  will  be  referred  to, 
henceforth  as  DNET-I. 

Additionally,  these  18  problems  were  chosen  since  they  are  gener¬ 
ally  available  to  other  researchers  and  the  problem  set  is  reasonably  small. 
The  specifications  of  these  18  problems  as  required  on  the  input  cards  to 
the  network  generator  in  [19]  are  given  in  Table  1.  Problems  1-5  are  100  x 
100  transportation  problems,  problems  6-10  are  I5O  x  150  transportation 
problems,  problems  11-15  are  200  x  200  assignment  problems,  and  problems 
16-18  are  1000  node  transshipment  problems.  Table  II  contains  the  solution 
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times  for  each  of  these  problems  for  each  of  the  codes. 

Comparing  the  advanced  start  dual  code  DNET-I  against  the 
previous  dual  code  DNET  indicates  that  our  initial  premise  (i.e.,  since 
the  majority  of  the  computational  time  spent  in  the  dual  procedure  is  de¬ 
voted  to  the  selection  of  an  in-coming  arc,  a  considerable  amount  of  time 
could  be  devoted  to  finding  the  starting  basis,  if  the  start  provided  a 
substantial  reduction  in  the  number  of  pivots  to  reach  optimality)  is  correct. 
In  particular,  the  results  show  that  the  difference  in  start  times  for  the 
two  codes  is  negligible  and  in  some  cases  (surprisingly)  the  advanced  start 
takes  less  time  to  find  the  starting  Dasis.  In  terms  of  pivot  reduction, 
the  advanced  start  is  quite  extraordinary.  For  example  at  least  a  1^  re¬ 
duction  is  achieved  on  every  problem  and  in  five  cases,  the  reduction  is 
more  than  'j&jo.  This  pivot  reduction  is  directly  reflected  in  the  total 
solution  time  since  DNE'T-I  strictly  dominated  DNET  in  each  case  and  the 
percent  reduction  in  total  solution  time  is  almost  identical  to  the  percent 
pivot  reduction.  Unfortunately,  this  substantial  reduction  is  not  enough 
to  make  the  dual  method  competitive  with  PNET-I  or  SUPERK  on  this  broad  class 
of  problems.  However  DNET-I  is  comparable  with  tbe  other  codes  Bennington, 
SHARE,  Boeing,  General  Motors,  and  Texas  Water  Development  Board.  Another 
noteworthy  feature  of  the  computational  results  is  that  PNET-I  strictly 
dominates  all  of  the  codes.  In  fact  PNET-I  is  roughly  twice  as  fast  as 
SUPERK,  the  fastest  known  out-of-kilter  code.  This  result  is  extremely 
important  since  it  completely  contradicts  the  folklore  that  the  out-of-kilter 
method  is  the  fastest  (in  terms  of  total  solution  time)  means  for  solving 
transportation  and  transshipment  problems. 

Motivated  by  Harris'  [19]  hypothesis  that  the  dual  simplex  method 
might  be  an  effective  solution  approach  for  solving  extremely  rectangular 
transportation  proolems  (i.e.,  problems  which  have  only  a  few  origins  and 
many  destinations),  we  decided  to  solve  some  rectangular  problems  using 
the  fastest  code  of  each  type.  In  order  to  make  our  research  results 
available  to  fellow  researchers,  we  used  NETGEN  [19]  to  generate  12  highly 
rectangular  transportatior  problems.  Table  III  contains  the  specifications 
for  these  problems  as  required  on  the  input  cards  for  NETGEN.  The  computa¬ 
tional  results  on  these  problems  are  given  in  Table  IV. 

Important  results  which  can  be  gleaned  from  Table  IV  are  that 
DNET-I  takes  less  time  than  SUPERK  in  all  but  one  of  these  problems;  how¬ 
ever,  DNET-I  again  assumes  its  secondary  role  to  PNET-I,  being  slower  on 
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all  of  the  problems.  Even  though  DNET-I  is  slower  than  PNET-I  on  these 
rectangular  problems  the  magnitude  of  the  differences  is  quite  small  corn, 
pared  wild,  the  fable  II  results.  Thus  Harris'  hypothesis  appears  to  have 
been  well  founded. 

Unlike  the  square  problems  of  Table  II  in  which  start  selection 
time  is  a  fairly  small  proportion  of  the  total  computation  time,  tectangular 
problems  cause  the  dual  to  spend  a  great  deal  of  time  picking  the  start. 

For  example,  the  problems  where  DNET-I  performed  best,  in  particular, 
problems  1,^,5, 7,  and  11,  the  start  (tine  constituted  9^$,  82$,  71$ , 

76$,  and  79^  respectively  of  the  total  computation  time.  This  appears  to 
reinforce  our  original  premise  that  quite  a  bit  of  time  could  be  devoted 
to  start  selection  if  that  start  reduces  the  numbev  of  pivots  being 
performed. 

More  research  is  needed  in  the  area  of  determining  precise  bound¬ 
aries  within  which  DNET-I  performs  best.  Initial  results  indicate  that  it 
performs  better  as  the  number  of  6*urces  decreases,  but  there  are  also  trends 
indicating  that  arc  density  plays  an  important  role.  Problems  8  and  9  show 
that  as  the  arc  density  increases  from  90$  to  90$  the  number  of  pivots  re¬ 
quired  to  reach  an  optimal  solution  nearly  doubles.  Yet  in  problems  3  and 
h,  an  increase  from  roughly  99$  to  roughly  75$  caused  nearly  a  30$  reduction 
in  the  number  of  pivots.  In  contrast,  PNET-I  is  not  as  sensitive  to  arc 
density  ranges,  since  the  number  of  pivots  is  proportional  to  the  number 
of  ares. 
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TABLES  IV 

Solution  Tines  (in  seconds) 


unht  -  i 


Problem 

Time 

#  Pivots 

1 

.292 

1 

2 

1.960 

58 

3 

10.249 

182 

4 

9.656 

122 

5 

18.996 

43 

6 

22.214 

88 

7 

21.866 

65 

8 

32.717 

261 

9 

56.926 

493 

10 

58.106 

313 

11 

43.735 

90 

12 

124.822 

758 

PNRT  - 

Time  # 

I 

Pivots 

SUPERK 

Time 

.079 

18 

.882 

.447 

84 

3.397 

2.414 

203 

13.330 

2.782 

139 

12.424 

3.927 

129 

35.287 

4.557 

128 

36.721 

2.734 

60 

37.547 

5.852 

259 

39.104 

10.314 

338 

43.726 

t' .286 

318 

67.459 

5.553 

114 

71.042 

18.024 

424 

DNR 

DNR  -  DID  NOT  RUN 
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